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Phenomenological Theory of a Scalar Electronic Order: 
Application to Skutterudite PrFe4Pi2 

Annamaria Kiss* and Yoshio Kuramoto^ 

Department of Physics, Tohoku University, Sendai, 980-8578 

By phenomenological Landau analysis, it is shown that a scalar order parameter with the point- 
group symmetry Fi^ explains most properties associated with the phase transition in PrFe4Pi2 at 6.5 
K. The scalar-order model reproduces magnetic and elastic properties in PrFe4Pi2 consistently such 
as (i) the anomaly of the magnetic susceptibility and elastic constant at the transition temperature, (ii) 
anisotropy of the magnetic susceptibility in the presence of uniaxial pressure, and (iii) the anomaly 
in the elastic constant in magnetic field. An Ehrenfest relation is derived which relates the anomaly 
of the magnetic susceptibility to that of the elastic constant at the transition. 

KEYWORDS: skutterudite, PrFe4Pi2, scalar order, magnetic susceptibility, elastic constant, 
uniaxial pressure 



1. Introduction 

Rare earth filled skutterudites RT4X12 attract continuous interest because of their intriguing and 
rich variety of phenomena depending on the rare earth R. Owing to intensive experimental and the- 
oretical efforts, some of these phenomena are now understood and explained, but there still remain 
unsolved problems of great significance. Among them are mysterious ordered phases in PrFe4Pi2/''^ 
metal-insulator transition and magnetic transitions under high pressure,^ and ferromagnetic ordered 
phase when a few percent of the Pr concentration is substituted by La.^ The purpose of this paper is 
to provide a coherent phenomenological picture for the fundamental order in PrFe4Pi2, in terms of a 
scalar electronic order. Other ordered phases in PrFe4Pi2 have different order parameters, and require 
separate analysis. 

PrFe4Pi2 undergoes a phase transition at Tq = 6.5K, which accompanies a sharp peak in the 
magnetic susceptibility.^ Staggered dipole moments are found in the ordered phase in the presence of 
a magnetic field. This observation indicates that the order parameter does not break the time-reversal 
symmetry.^ A modulation with the wave vector Q = (1, 0, 0) was found below Tq by X-ray diffraction 
experiments,^ which is attributed to the existence of staggered local electronic states of the Pr ions. 
Early NMR^ and elastic measurements*^ were interpreted in terms of antiferro ordering of r3-type 
quadrupole moments, which picture was widely accepted for a long time. However, the isotropy of 
the magnetic susceptibility for different field directions even in the ordered phase cannot be explained 
by quadrupolar order. Furthermore, staggered dipoles perpendicular to the field direction are found 
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neither in neutron diffraction^'^ nor NMR.'" Recently, careful analysis of the NMR results pointed 
out that the local symmetry at the Pr sites is preserved in the ordered phase. ' ^ Furthermore, the 
continuous field-angle dependence of the transition temperature gives also evidence for the exclusion 
of r3 quadrupolar order. 

In a previous paper, we proposed that the order parameter in PrFe4Pi2 is a staggered electronic 
order which does not break the local symmetry around each Pr site. We call this order with the 
Tig symmetry a scalar order. It is found in ref. 13 that the scalar order scenario can explain naturally 
the isotropic magnetic susceptibility in the ordered phase, the field angle dependence of the transi- 
tion temperature and magnetization, and also the splitting pattern of the ^^P NMR spectra. The main 
purpose of this paper is to provide a coherent description of available magneto-elastic properties of 
PrFe4Pi2 within the scalar order scenario. We can explain not only elastic anomalies in PrFe4Pi2 near 
the phase transition, but also the huge anisotropy of the magnetic susceptibility under uniaxial pres- 
sure. Our model assumes phenomenological couplings of the scalar, dipolar, and quadrupolar degrees 
of freedom. 

This paper is organized as follows. In section 2 we introduce the free energy expansion around 
the transition temperature and clarify its constituents. Section 3 discusses the properties of magnetic 
susceptibility, and its anisotropy in the presence of uniaxial pressure is obtained. In addition, particular 
attention is paid to the emergent anisotropy in magnetic field. In section 4 we discuss elastic properties 
of the scalar order within the phenomenological framework. Section 5 studies Ehrenfest relations 
relating to the anomalies of magnetic and elastic quantities. The last section is devoted to the summary 
of the paper. 

2. General framework of phenomenology 

2.1 Landau expansion of the free energy 

A scalar order can be two different kinds: one which breaks the time reversal symmetry (Fij,), and 
the other which does not (Fi^). In this paper we assume a two-sublattice scalar order with the ordering 
vector Q = (1, 0, 0), where the order parameter has the T\g symmetry. 

In the Landau theory, one expands the free energy in terms of the set of electronic order parameters 
Mi*,, which are taken to be real. Up to fourth-order, we write 



where we have introduced the quantity a,- = a,(r - T,). 

The constants a,-, bi are positive, while gij and Cij can have either sign. T, is a hypothetical transi- 
tion temperature without coupling to other order parameters. The actual transition occurs at Tq corre- 
sponding to the scalar order, which gives the largest of all T, . For other component ^„ we neglect the 
corresponding bi in most cases. These parameterizations have a merit that each coefficient is regarded 
as a constant as long as the temperature is close to Tq, and the external perturbations are small. 

As explicit constituents of we include the scalar order parameter ij/Q, the homogeneous mag- 
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netization M, and the Fs-type homogeneous quadrupoles. Furthermore, we also include the lattice 
strain components Sxx, £>•>■» ^zz^ which have a bilinear coupling with quadrupole moments Qij. The 
second-order couplings with / = Q, ipQ can be neglected because the background elastic constant 
C^°^ is large enough. The external magnetic field is assumed with general direction, and the magneti- 
zation components are included up to fourth order. We write the free energy consistent with the cubic 
symmetry as 



r(i/fQ,Q,M,s) = 



2 



+ 



8mq 



[iMl - Ml - Ml) Qu + ^ (Ml - m;) 



(2) 
(3) 



where = a^(T - To), oq = agiT - Tq), um = aM{T - Tf) and M^ = Ml + M^ + M^. For the strain 
components we introduce the notations £„ = (1/ ^^6)i2£^y - e„ - £„,), = (1/ ^|2){Exx - %>) and 
= (1/ V3)(£xx+ejT+ezz), and cfl and C^^^^ are the elastic constants in the background. Qi, and 2v are 
the F3 quadrupole moments defined as e„ = = (1/ ^l6){2Jl -JI- J^), Q^ = 0\ = {\I ^){JI - Jy), 
and we introduced the notation 2^ = 2^ + 2v- 

As external fields, we consider only two cases: magnetic field and uniaxial stress. We will consider 
the uniaxial stress direction (t||(001), so that we have cr^^ + and cTxx - CTyy = 0. This stress can be 



decomposed into the sum of an isotropic and a traceless part as 



cr 



yy 





cr 



1 



-1 
-1 
2 



(4) 



Here, the first part (1, 1, 1) has Fi symmetry and it is the isotropic hydrostatic pressure component. The 
second part (-1,-1,2) has F3 symmetry, and it causes the lowering of the original cubic symmetry 
to tetragonal one. For the uniaxial stress components we introduce the notations cr„ = (1/ ^f6)(2(r^^ - 
o-xx - o-yy), (Tv = (1/ ^XcTxx " (Tyy) and OS - (1/ ^/3){(rxx + (ryy + a-^), where we have the constraint 
(Tu = V2o"i. and cTy = 0. 

The equilibrium condition leads to the coupled equations: 



dr 

^ - a!M,kMk, 

oMk 

o"m = ^ = C^Su + BQu, 

dr 

0^ — = C3Sy + BQ„ 

OEv 



dr -> 

O", = ^ = CqEs + gtl,s^Q + gMsM^ + gQsQ^, 



&r 1 / , n 

= ^ = ccqQu + Se„ + —guQ (2Mj - Ml 

dr 1 / 9 7\ 

= ^ = aQQv + Bsy + -j^gMQ [Ml - M;) , 



Ml), 
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= -r— = a^^ii, (5) 

where k - x,y, z, and a,- indicates a renormalized inverse susceptibility. They are given together with 
C3 and Co by 

aM,;c = Q^M + 2gMs£s + C^M^lfQ + CmqQ^ + gMQ-^ ( VSQv - Qu) , (6) 

2 

Q^Mj -aM + 2gMsSs + C^M^i^ + CmqQ^ + gMQ-^ (" VSQv - (2m) > (7) 

2,4 

aM,z -0£M + 2gMESs + C^M^i^ + cmqQ + gMg— (8) 

V6 

Q;j2 = + 2g2e£i + c^q^/'q + cmqJ^^, (9) 
a^ = a^ + Ig^sEs + b^ij/^ + c^qQ^ + c^mJ^^, (10) 
C3 = C-C(?, (11) 
Co^cfi^ + ici^. (12) 
2.2 Phase boundary 

The transition temperature = To + ATq is obtained by the condition - in eq.(lO). Let us 
first consider the case under magnetic field, but with zero stress. We assume that the magnetic field is 
parallel to the direction (001), which gives -My-Q, and we write = M. We note that the final 
result does not depend on the field direction up to second order in the field. The equilibrium condition 
is given by 

H = aMM, (13) 
= Cos. + gMsM'^, (14) 
= a^il/Q - (a^ + c^mM^ + 2g^sSs)^Q- (15) 

The condition - gives the temperature-magnetic field phase boundary. By eUminating M and 
in eq.(15) with the aid of eqs.(13) and (14), we obtain 

We will neglect the terms where the background elastic constants Co or C3 appear in the denominator 
since these terms are small. The exception arises when these terms are multiplied by the uniaxial 
pressure, because the product (l/Co)cr is not small. Thus, we obtain from expression (16) 

a^(%{TQ - TfY 

From the measured temperature-magnetic field phase boundary we obtain tHi ~ 0.34[K-T~-^].^ There- 
fore, we estimate Cfi,Ml[a,f,cP^{TQ - Tpf-^ ^ 0.34[K-T~-^]. 

The phase boundary under uniaxial stress cr ^ can be obtained in a similar manner. The ex- 
perimental temperature-uniaxial pressure phase boundary is linear for low values of cr. Hence we 
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Fig. 1 . Temperature-magnetic field (left) and temperature-uniaxial pressure (right) phase boundaries close to 
the transition temperature Tq = 6.5K. For the parameter values, see the text. Black boxes represent the 
measured result taken from refs. 1 and 14. 



concentrate on the linear regime, where only the scalar combination cr^ = 1 / ^{(Txx + (Tyy + cr^) 
enters. Without magnetic field, we obtain 

Tc^To —o-s = To- to-io-s. (18) 

a^Co 

The estimate ta-i/ ~ 1.5[K-GPa~^]^'^ gives the restriction for the combination: g,j,E/(a,j,Co) ~ 
1.3[K-GPa~^]. The fitting result for the T-H and T-cr phase boundaries are shown in Fig. 1. 

For later convenience in discussing Ehrenfest relations, we also go on to the second order regime, 
where two different contributions enter: cr^ - l/3((rxx + o-yy + cr^ and (tI + crl - 2/3(o-^^ + cr^ + 
(T^ - (TxxO'yy - (TxxO-zz ~ o'yyf^zz)- Since we neglected the second-order couphng c^g in the free energy 

ryy 

we obtain 

1 



expansion (2), the term with ((Txx + o-yy + cr^^)-^ cannot be calculated. For the other second order term 



[c,a - i^l + crl) . -t^, (al + crl) . (19) 



[S2 - ag(ro - TQ)f- 
3. Magnetic susceptibUity 

3.1 No uniaxial stress 

In the case of zero uniaxial stress, the magnetic susceptibility is obtained from eq.(2) by setting 
(Tu - (Ts - O va the equihbrium condition (5). For T > Tq, we further set if/Q - 0, and obtain the 
Curie- Weiss law: 

X-.'-q^2-^m{T-T,), (20) 

where Tp is the Weiss temperature, and we use the superscript + to indicate T > Tq. In the ordered 
phase close to the transition temperature, we obtain 



c,/fMai/f(T - To) 2a^(T - To)gi/,E I c^Mg^e 



= auiT-Tp) -^—TT^ Tl — [s^e r > (21) 
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where we set t/^g and e., using the equilibrium conditions &F/dif/Q = and Idss = 0. Neglecting 
the last term in eq.(21), we obtain the inverse magnetic susceptibility as 

;^:^-aM(r-r.)- ^^""^f-^"\ (22) 

There is a peak in the magnetic susceptibility dlT -Tq if dxZ^ /dTlrg < is satisfied. The temperature 
derivative is calculated as 



= aM 7 — , (23) 



dT 

which gives the condition a^c^M > ciMb,p for occurrence of the peak 'mx{T) 31 T = Tq. Expressions 
(20) and (22) for the magnetic susceptibility are used to fit the measured result. We set Tp - 3.5K 
for the Weiss temperature, which is found in the experiments.^ The value for the parameter qm is 
obtained from the fit in the paramagnetic phase as um - 9.6 • 10^[Pa-yUg^-K"^]. Additionally, the 
combination a^fjC^Mlbip is obtained from the fit of the susceptibihty in the ordered phase. We can 
choose the value of b^ and then our previous estimation for the parameter tm from the measured 
temperature-magnetic field phase boundary gives the values of and c^m- The results are given by 

= \0'^[Pa], ^ 1.95 • 10'^[Pa-K-i], c^m = 1-37 • 10'^[Pa -jUg^]. The fitting result obtained is shown 
in the left part of Fig. 2. 
3.2 Finite uniaxial stress 

Intriguing experimental results are obtained for the magnetic susceptibility in the presence of 
uniaxial stress. ^^'^^ Namely, the magnetic susceptibility shows large enhancement for the uniaxial 
pressure applied parallel to the magnetic field direction {H\\o-), while it shows only slight decrease 
when the pressure is applied perpendicular to the field direction {H ± cr). 

Now we discuss the properties of the magnetic susceptibility around the transition temperature Tq. 
The direction of the uniaxial stress is taken as (t||(001), and we consider two different directions of the 
magnetic field, namely //||(001) and //||(100). These two cases can be obtained from the free energy 
given by eq.(2) by taking M^i^O and + 0, respectively. For small values of the uniaxial stress, it is 
enough to consider only the Unear term in cr. Thus, we obtain the susceptibiUties x\\ for //||(001) and 
X^ioTH\\{\m) as 

-1 4 2 

X\\ -aM+ -^EMqQu + '^gMeSs + C^M^Q^ (24) 

X~i} - am- -^SmqQu + '^gMe^s + c^Mt/fq- (25) 

Among the terms appearing in the expression of the magnetic susceptibility, gMs^s and SmqQu con- 
tain the uniaxial stress. We can find in eqs.(24) and (25) that the former term is isotropic, while the 
latter term is anisotropic and has different sign for the two magnetic field directions. Therefore, the 
cancellation of the isotropic and anisotropic terms can occur for the case H ± cr, which reproduces 
the experimental behavior. 

Now we set i/fQ, Ss and Qu from the equilibrium conditions dTldif/Q - 0, dTldss - os - 
Vl/3cr and dT^/dsu - (Tu - ^/2/3o■. Keeping only the leading term as we did also previously in the 
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Fig. 2. Left; Magnetic susceptibility around the transition temperature Tq. Black boxes represent the measured 
result taken from ref. 1. Right: Uniaxial pressure dependence of the susceptibility peak at T = Tc for the 
magnetic field directions (001) and (100). Black boxes represent the measured result taken from ref. 14, 
where we converted the magnetic susceptibility from unit of [emu-mol~^] to [/t^-ion"^ T"^]. 



susceptibility calculation with no uniaxial stress, we obtain the magnetic susceptibilities x\\ and;\fj_ in 
the ordered phase as 

-1 -1 4 / BgMQ \ 2 IgMs] 

-1 -1,2/ BgMQ \ 2 (gMs\ 

where is given by eq.(22), and 

c^Qa^jT - To) 

aq-QQ . (28) 

The susceptibiUties in the paramagnetic phase can be obtained by taking oq - oq in eqs.(26) and 
(27). In experiment, the uniaxial pressure dependence of the magnetic susceptibihty at the transition 
temperature^^ is Unear for small values of cr. We make series expansion in terms of cr in the expressions 
(26) and (27), and obtain the susceptibiUties x\\ and;^fj_ at the transition temperature Tc as 

where we introduced the notation 

y.J ^ (31) 

^B^-C^aQ{Tc-TQ) 
and Tc means the transition temperature in the presence of uniaxial pressure given by eq.(18). 

If the condition gMs/Co ~ 7 is satisfied accidentally, the term with uniaxial pressure cr disap- 
pears in the expression of x±(Tc)- We interpret the experimental situation in this way, and obtain the 
susceptibilities as 

XwiTc) = — „ . + r • 0-, X±iTc) = — rr (32) 

aM(Tc - Tf) aM(Tc - Tp) 
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Fig. 3. Field angle dependence of function f{6, <p) with (f> = n/4 (left) and 6 = 7t/2 (dght), which means the 
planes [1,-1,0] and [0,0,1], respectively. The parameter values are chosen ast = 1/2 and = -1. 



where 

a^iTc - Tp)^ fi2 _ _ J^^c^ a^iTc - Tp)^ Co " 
We find that with the choice of F = 2.85 • 10^[j[i|-GPa~^] the measured susceptibilities xwiTc) and 
Xi-iTc) can be fitted at the transition temperature. The result is shown in the right part of Fig. 2. 
Different anisotropic behavior is reaUzed depending on the parameters gMe and gMQ- 

From the expression (2) we can calculate the magnetic susceptibihty for a general direction of the 
magnetic field with respect to the uniaxial pressure cr||(001). We obtain the magnetic susceptibility for 
the magnetic field direction {hx, hy, h^) at the transition temperature as 



aM(Tc - Tp) 



ai (hi + hy+ hfj + a2^ (hi + hj - 2hj^ 



(T, (34) 



where 



^ 2 I gMs ^ ^ 4 1 BgMQ 

V3 4(r, - Tp)^ Co ' 3 al{T, - Tp)^ {B^ - C,aQ(T, - Tq)) ' 

We require that in the case of // J_ cr with h^ = I, hy = h^ - 0, the uniaxial pressure term disappears 
in the magnetic susceptibility. Using this condition ai + l/2a2 = 0, we obtain 

XiTc) = ^ , + ^aihjo- = ^ + (3ai cos^ o) a. (36) 

3.3 Nonlinear effects of magnetic field 

In order to discuss the anisotropy under magnetic field, we have to include the fourth order terms 
of magnetization in 'F{if/Q,M). It is then more convenient to perform the Legendre transformation to 
magnetic field 7?, to obtain the free energy expansion as 

n>I^Q, H) = ro{H) + \a^[T - TciHM^Q + ^^^1^^, (37) 

where 

Tc(H) = To - tH2H^ + (tA + t4ah4)H'^ (38) 
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is the transition temperature in the presence of magnetic field. Here, tni, k and t4a are coefficients, 
and h4 is the fourth order cubic invariant h4 - h'^ + hy + hf - 3/5, which is common in both cubic Oh 
and tetrahedral J/, symmetries. The reference part Toili) has the field dependence similar to Tdfl). 

The second order term in the expression of the transition temperature given by eq.(38) is isotropic, 
while the fourth order term carries anisotropy expressed by the function f{d, (p) = t + tah^, where the 
field angle Q and (p are defined as {h^, hy, h^) - (sinf^cos^, sinf?sin0, cos0). We note that the same func- 
tion f{Q, (p) determines the magnetization anisotropy as well for small magnetic fields. Interestingly, 
the anisotropy given by the function f{6, (p) is independent of the parameter values of the free energy 
expansion, i.e., the microscopic details of the scalar order. Figure 3 shows the field angle dependence 
of the function f{G, (p) in the [1,-1,0] and [0, 0, 1] planes. We find that for the three principal axes of 
the magnetic field, namely for the field directions (001), (110) and (111), the following ratio 

mn-m^^ - 4 (39) 

/[110]-/[111] ^ ' 

is held, where we used that ;z4[001] = 2/5, /j4[110] = -1/10 and ;j4[lll] = -4/15. Thus, the ratio 
(39) is satisfied for the transition temperature and also for the magnetization as long as the magnetic 
field is small. The function f{6, <p) gives an excellent fit for both the observed transition temperature 
and magnetization in the case of PrFe4Pi2 with proper choice for the parameters t and ta, which gives 
strong evidence for scalar order in this compound. 

4. Elastic constant 

4.1 Case of zero magnetic field 

The free energy expansion in zero magnetic field can be expressed in terms of the strain compo- 
nents as 

ns) = ^CijkiSijSki, (40) 

where the elastic constant Cij^i can be obtained as Cij^i = d^T I dsijdeki. The connection between the 
strain and stress components cr,-y is given by cr;^ = ^ Cijusu- 

We discuss the properties of the elastic constants Cn - Cu which is related to the Tt, quadrupoles, 
and Cii + 2C\2 which is related to the scalar order parameter. They can be obtained from the free 
energy expansion (2) as 

Cn-Cn^^^^{CzEu + BQu), Cn +2Ci2 ^ ^ ^ (CqSs + 2g^s^Q) . (41) 

Let us consider first the elastic constant Cn - C12. The condition d'F/dQu - gives that Qu - 
-{BlaQ)Su, where chq contains and e,. We set i/^q and Ss using the equilibrium conditions 
dT'/diJ/Q - and d'F/dsg - 0, and obtain 



Cn - C\2 - C3- 



b^Co-2g^ \ b^ 
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Fig. 4. Left; Elastic constant given by the expression (43) around the transition temperature To in the case 
of zero magnetic field. Right: Elastic constant in magnetic fields H = 0, 0.5 and IT given by the series 
expansion (52). Black boxes represent the measured result taken from ref. 16. The parameter values are 
given in the text. 



Neglecting the last term of the denominator in eq.(42), we obtain 

Cii - Ci2 = C3 — = cf) - Cf^ t;—- (43) 

For the elastic constant Cn + 2Ci2 we obtain 

C„ + 2Cn = Co-^4^= Cf^ + 2dS ' (44) 
In the disordered phase (^^^q = 0), the elastic constant Cn - C12 is modified as 

B2 B^xq 



SiCn - Cn) = - 



where we introduced 5Cij - Cjj - C^^\ ^^xq is the quadrupolar susceptibility. At high temperatures 
it has the property (5(Cii - Cn) ~ ^/T, which is analogous to the Curie law. 

In the ordered phase close to the transition temperature we obtain the change 6Cij of elastic con- 
stant as 

6(Cn-Cn)- -z; — zttt t:^ — (45) 

If the condition d (Cn - C12) /STItq < is satisfied, the elastic constant Cn - C12 has a negative peak 
at T - To due to the scalar order. The temperature derivative of Cn - C12 atT - Tq is calculated as 



5(Cn-Ci2) 



dT 



^ d6iCn-Cn) 
To dT 



To b^ayjo - Tq)^ 
which gives the condition a^c^g > agb^ for the negative peak. 

We take the value for C3 from experiment, which gives Cfl - = VOGPa at r = 4.2K,^ and 
choose Tg - 0.4K. Then we obtain the values for the combinations B'^ lag and B'^/c^g from the fit of 
the experimental data^^ as B^/ag - 10[GPa-K], B^/c^g = 2.2[GPa]. The calculated elastic constant 
Cn - C12 around the transition temperature Tq can be seen in the left part of Fig. 4. 
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Around the phase transition we make series expansion of expressions (43) and (44) in terms of 
(T - To), and obtain the change dCtj of elastic constants as 



S(Cn - Cn) 



6(Cn + 2Cn) 



aQ{T - Tq) 

,2 



1 + 



a^jjCiiiQ T — Tq 



b^UQ T -T, 



Q 



2gf 



(47) 



(48) 



In the present mean field theory 5(Cii - C12) is continuous at T = Tq, and only its temperature 
derivative has a discontinuity. On the other hand, the elastic constant Cn+lCu is discontinuous at T = 
To, namely A(Cii + 2C12) = ^(Cn +2Cu)Tg -SiCn +2Ci2)t^ * 0. and the jump A(Cii + 2C12) 
is always negative in the mean field theory. 

Using the values for the parameters obtained, the change of the elastic constant Cn - Cn around 
the transition is estimated as 



AC 
AT 



diCn-Cu) 



dT 



To 



2.3GPaK-\ 



(49) 



^uq) S2 b^{To - Tq)^ 

which has the same order of magnitude as the inverse of parameter ta-i- The latter describes the 
temperature-uniaxial pressure phase boundary as to-i - ATc/Acr (see eq.(18)). We also estimate the 
relative change in the elastic constant Cn - C12 around the transition temperature as 



1 AC 
CAT 



1 



^(0)_^(0)' 
'-'11 '-'12 



^(Cll-Ci2) 

dT 



3.2 - 10"^ K-\ 



(50) 



T=To 



which means a few percent change of the elastic constant around Tq compared to the background value 



^1? ~ ^12' accordance with the experimental situation. 
4.2 Case of non-zero magnetic field 

Now we consider the effect of non-zero external magnetic field with direction (001) on the elastic 
constant Cn - C12. We use the relation 



-(0) ; 



Cn - C12 = — = ^ (C3£u + BQu) . 



(51) 



The condition d'F/dQu - gives the quadrupolar moment as Qu - -(B/aQ)su, where oq now 
contains also the magnetization M^ = M besides iJ/q and e^. We set ^ and Ss using the equilibrium 
conditions dTldif/Q - 0, dTldss - 0, and H - cumM. We make series expansion in terms of T - Tq 
in the vicinity of the transition temperature Tq, and also in terms of the magnetic field since it appears 
in the denominator. We obtain the elastic constant as 



Cn - C12 



S2 

C3 + 



B' 



CMQ + T 



//^-^^(r-r.) 



+ 



2B^ 



3 aM ^ 

CipM _^ C^Q \ I C^mC^Q \ ^ 2glfQ / c^M _^ 2 C^jjQ \ 

am ag j b^ j aM [cum 3 / 



b^ 

2gl 



Q "f 

2 



(T - Tc)H\ 
(52) 
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Fig. 5. Phase boundary in the T -X plane around the transition temperature To . 



where Tc means the transition temperature in the presence of external magnetic field given by eq.(17). 
In the paramagnetic phase the elastic constant reduces to 



Cii - Ci2 - C3- 



■>2 



4 8mq 

CMQ + T - 



3 UMiT - Tp) 



h\ (53) 



aQ(T-TQ) aliJ-TQfaliT-Tff 
which shows a softening of the elastic constant in magnetic field provided the coupling constant cmq 
has a negative sign. Furthermore, using a proper choice for cmq, the coefficient of term {T - Tc)H^ in 
the ordered phase can be negative. Then the sharpening of the anomaly occurs at Tc in magnetic field 
since the absolute value of the derivative d (Ci 1 - C12) /dTlj^ increases. Right part of Fig.4 shows the 
behavior of the elastic constant in magnetic field with the parameter value B^Jcmq - -0.5[GPa-/ig]. 

5. Ehrenfest relations 

Let us regard a continuous phase boundary on the T - X plane, where T is the temperature and 
X is an arbitrary external quantity such as uniaxial stress, hydrostatic pressure or magnetic field (see 
Fig.5). The vicinity to the phase transition Tq is considered, where the quantity X is small. We assume 
that the phase transition is of second order around Tq. The free energy is continuous along the phase 
boundary, i.e. T{Tq + dT, dX)i - T{Tq + dT, dX)2, where 1 and 2 labels the ordered and paramagnetic 
phases, respectively. Furthermore, in the case of a second order transition the first derivatives of the 
free energy dTjdT and dTjdX are also continuous along the phase boundary. Two difi'erent cases are 
considered: (i) Unear phase boundary given by Tc - Tq- txiX and (ii) quadratic phase boundary given 
by Tc - To - txiX^. Using the continuity of the free energy and its first derivatives along the phase 
boundary, we can obtain relations between the anomalies of difi'erent thermodynamic quantities at the 
transition temperature Tq. These relations are called Ehrenfest relations. 
5.1 Linear phase boundary 

In the case where the phase boundary has a Unear dependence such as Tc - Tq- txiX, we obtain 
the following Ehrenfest relations at the transition temperature Tq 
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, AC, d'^T 

where is the specific heat defined as = -Td^T jdT^. 

With X being the uniaxial pressure {X - cr), the T-cr phase boundary was obtained in eq.(18) for 
small values of the uniaxial pressure as Tc - Tq - ta-\crs, where cts - 1/ ^((Txx + o-yy + cr^ is the 
hydrostatic pressure. The quantity d^TjdcrsdT = or in the Ehrenfest relations (54) is called thermal 
expansion coefficient a. Using the equiUbrium conditions d'F/diJ/Q - and dTldss - <Ts where the 
free energy T is given by eq.(2), the discontinuity in a at the transition is calculated as 

b^Co - 2g^^ 

To obtain the anomaly of the specific heat, we calculate the temperature dependence of the free energy 
(2) in the absence of external fields. Using the equilibrium conditions dTldif/Q - and dTldsg - 0, 
we obtain 



Aa = 7"^"^ ^ . (56) 



1 Coa^(r - ToY 
4 b^Co-2gl 



which gives the jump of the specific heat atT -Tq as 

To ~2b^^ (fe^Co - 24) " 2 b^Co - 24 • 
Using eqs.(56) and (58) together with the expression of ta-i given in eq.(18), we find 

-f^i^-Aa = 0, (59) 

which is the special case of (54) with X = cTg. 

Since the specific heat shows a jump at the transition in PrFe4Pi2/ we expect discontinuity also of 
the thermal expansion coefficient. However, there is no observed discontinuity of the thermal expan- 
sion coefficient in ref. 17. Let us compare with the case of URu2Si2 where experiment with hydrostatic 
pressure has been performed. ''^ In the case of URu2Si2, both the specific heat and thermal expansion 
coefficient show a jump at the hidden order phase transition, and the Ehrenfest relation (54) is satisfied 
within the experimental error. The observed discontinuity of the thermal expansion is about an order 
of magnitude smaller than the error bar in the experiment for PrFe4Pi2.'' Both URu2Si2 and PrFe4Pi2 
have almost the same value in the combination t^i AC JTq in eq.(54). Therefore, in PrFe4Pi2, we expect 
similar magnitude of discontinuity of the thermal expansion coefficient at the transition temperature, 
which may be observed with improved resolution. 

Let us discuss the Ehrenfest relation (55). The compliance +2512 given by -d^f/dX^ with 
X = cr,; in eq.(55) is related to the elastic constant Cn + 2Ci2 as 5 n + 25 12 = (Cn + 2Ci2)~^ From 
expression (44) we obtain the discontinuity of the elastic constant Cn + 2Ci2 as 

A (Cii + ICn)-' = ^^^-TT' (60) 

Co(b^Co-2gl^) 

where A(Cii + 2Ci2)-^ = (Cn + 2Cii)-\ - (Cn + 2Ci2);l 
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tli^-MCn+2Cnr' ^0, (61) 



Using eqs.(58) and (60) together with the expression of ta-i given in eq.(18), we find 

AC, 
To 

which is the special case of (55) with X - <Ts, and relates the jump of the elastic constant to that of the 
specific heat. Since the jump of the latter is finite at the transition temperature in the mean field theory, 
the inverse elastic constant (Cn + ICnT^ should also have a finite jump. However, this jump has so 
far not been found in the experiment in PrFe4Pi2.^ We beUeve that such a jump is indeed present, but 
is too small to be observed. 
5.2 Quadratic phase boundary 

When the phase boundary has the form Tc - TQ-txiX^, we obtain the following Ehrenfest relation 
at r = To in the lowest order 

AC, 1 d^T 

?x2^ + 2^Q^2gi - ^^^^ 

Let us consider the case where X is the external magnetic field {X = H). The temperature-magnetic 
field phase boundary is obtained as quadratic in the expression (16), namely Tc = Tq-IhiTI^- From the 
expression (21) we obtain the discontinuity of the temperature derivative of the magnetic susceptibility 
at Tq as 



(63) 



dT UMiTo - Tpf \ b^Co - 2gl^ 

Using the expression of parameter tH2 given by eq.(16) together with the jump of the specific heat 
given in eq.(58), we find the relation 

AC, 1 dy 

which is the Ehrenfest relation (62) with X - H. 

Similarly to the magnetic susceptibiUty, the elastic constant Cn - C12 also possesses anomaly at 
the transition temperature. Using the expression (42), we find that the discontinuity of the temperature 
derivative of Cn - Cn is related to the specific heat jump through the phase boundary parameter tcri 
given in eq.(19) as 

Ac,_i a(c..-CH)- 

Tq 2 dT 

which is the Ehrenfest relation (62) withX = cTu- Combining the expressions (64) and (65), we derive 
the Ehrenfest relation 

-La^<^"-^"'" = J-a|^. (66) 

ta2 OT tH2 OT 

which relate the anomaly of the elastic constant Cn - C12 to that of the magnetic susceptibiUty at the 
transition. We note that parameter t^-i gives the second-order term of the phase boundary relating only 
to the traceless part of the stress tensor (see eq.(19)). Therefore, it is difficult to extract tffi from the 
experimentally observed phase boundary. 

Finally, we summarize in Table I the roles of the parameters appearing in the free energy expansion 
(2) together with the restrictions for the parameters that we found in the Landau analysis. 
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Table I. Summary of the parameters appearing in the free energy (2). The restrictions for these parameters 
obtained from our analysis are also included. The parameter y is defined by eq.(31), and (t is the Weiss 
temperature 0* - 3.5/r.' 



Parameters 


Relevant quantities 


Restrictions 




X+iT) 


um >0,Tf = S* 




X-(T) 
T-H phase boundary 


C^M > 




T-(T phase boundary 


> 




(Cii-Ci2)(// = o,r>ro) 


flg > 




(c„-Ci2)(// = o,r<ro) 


ciil/C^Q > agb^ 


Smq, 8ms 


X\\{o-),xAcr) 
X±^X-iT = T,) 


EMe/Cq V y 


cmq 


{Cn-Cn)(H) 


Cmq < 



6. Summary 

In this paper we have shown that the main properties of PrFe4Pi2 can be explained consistently 
within the scalar order scenario. We have discussed the effect of a scalar order on magnetic and elastic 
properties, and appUed the results to the case of PrFe4Pi2. It is found that the anomaly of the magnetic 
susceptibility and the elastic constant at the transition temperature can be reproduced consistently 
together with the phase boundaries against to temperature-magnetic field and temperature-uniaxial 
pressure. Furthermore, we obtain the enhancement of the magnetic susceptibiUty for the uniaxial pres- 
sure applied parallel to the direction of the magnetic field, and the almost constant behavior of the 
susceptibility when the uniaxial pressure is appUed perpendicular to the magnetic field. In our view 
this behavior is caused by the competition of the isotropic and anisotropic stress dependent terms. The 
magnetic susceptibiUty is calculated also for arbitrary direction of the magnetic field with the uniaxial 
pressure fixed along (001). We have studied Ehrenfest-type relations which connect the anomaUes of 
different kinds of thermodynamic quantities at the transition temperature such as specific heat, mag- 
netic susceptibility or elastic constant. The phenomenological analysis leaves several questions open 
about PrFe4Pi2 such as the microscopic nature of the scalar order parameter, appUcabiUty of localized 
picture for the 4/ electrons or the ground state at low temperatures. Further studies are necessary to 
clarify the microscopic origin of the scalar order parameter in PrFe4Pi2. 
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